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We compare the decay of turbulence in superfluid "^He produced by a moving grid to the decay of 
turbulence created by either impulsive spin-down to rest or by intense ion injection. In all cases the 
vortex line density C decays at late time t a.s £ oc At temperatures above 0.8 K, all methods 

result in the same rate of decay. Below 0.8 K, the spin-down turbulence maintains initial rotation 
and decays slower than grid turbulence and ion-jet turbulence. This may be due to a decoupling of 
the large-scale superfluid flow from the normal component at low temperatures, which changes its 
effective boundary condition from no-slip to slip. 

PACS numbers: 47.80.Jk, 67.25.dk, 47.27.-i 


Turbulence is a common state of flow in classical flu¬ 
ids, with great importance from atmospheric systems to 
aircraft design. So far, satisfactory understanding is only 
achieved for homogeneous and isotropic turbulence (HIT) 
mm- HIT can be approximately obtained in the wake 
of a flow past a grid although it might still be 

strongly modified by the container geometry [HE!. Grid 
turbulence in superfluid "‘He was obtained [ZHS], but not 
at temperatures below 1K due to technical difficulties. 
Yet, the low-temperature regime enjoys a special inter¬ 
est, as the thermal excitations (the normal component) 
are essentially absent. Turbulence of the superfluid is 
made of a chaotic motion of tangled topological defects 
of the superfluid order parameter field - quantized vor¬ 
tices - each carrying the same circulation equal to the 
ratio of the Planck’s constant to the mass of a "‘He atom: 
K = hm~^. It is called Quantum Turbulence (QT), as it 
is essentially a macroscopic quantum phenomenon. QT 
decays even at the lowest temperatures, but the mech¬ 
anisms of dissipation in superfluid "^He - thought to be 
the radiation of phonons by Kelvin waves (perturbations 
of vortex lines) with wavelength ~ 10“^ cm (TUI and also 
of small ballistic vortex loops that can carry energy away 
[nHH - only operate at very small length scales. Ex¬ 
isting theories [MU of QT decay are applicable to ho¬ 
mogeneous isotropic QT (HIQT), for which only sparse 
experimental data are available in the interesting ultra- 
low temperature limit. 

In this Letter we report the best-yet realization of 
HIQT in the T —)■ 0 limit. We measure the free decay of 
grid turbulence and compare the results with both theory 
and experiments using other methods, thereby gaining 
valuable insights into the underlying processes. 

When QT is generated by large-scale flow, on length 
scales much greater than the mean intervortex distance 
Iq — where £ is the length of vortex lines per 

unit volume, then the energy is predominantly contained 
in flow at the largest length scales £q. In this case 
QT is called quasiclassical |22J [23], as quantization of 


vorticity becomes unimportant, and the coarse-grained 
velocity field is expected to obey the Euler equation. It 
is believed that this energy cascades towards the smaller 
length scales via a classical hydrodynamic cascade, fol¬ 
lowed, at length scales < £q, by a ‘quantum cascade’ 
that involves reconnections and Kelvin waves on discrete 
vortex lines. Existing theories [m dg dlHSO] of these 
processes in HIQT all assume that the dominant contri¬ 
bution to L is at quantum mesoscales ^ iq, but they 
differ in detail. For self-similar flows, assuming that the 
rate of dissipation of flow energy per unit mass, £, only 
depends on C and k, dimensional considerations demand 

£ = ( 1 ) 

Here, the ‘non-dimensional effective kinematic viscosity’ 
C ^ 1 (the more conventional ‘effective kinematic vis¬ 
cosity’ is v' = (k) [2g [23]. At medium temperatures 
1.0 < T < 1.6 K, it reflects the dissipation through the in¬ 
teraction of vortices with thermal excitations (expressed 
through the ‘mutual friction parameter’ a{T)), while in 
the limit T —>■ 0 (T < 0.5 K) it characterizes the effi¬ 
ciency of the tangle of vortex lines in maintaining the 
energy cascade down to the dissipative length scale. As 
there is no microscopic derivation of Q, it remains un¬ 
clear whether the value of ( is the same for HIQT of any 
spectrum, or depends on the type of flow. For instance, 
C = 0.08 was measured [25] for Vinen (‘ultraquantum’, 
i. e. without flow at classical length scales > iq) QT at 
T —)• 0, while the analysis of the decay of QT generated 
by spin-down at T —)■ 0 apparently revealed C « 0.003 
[26] . The latter was heralded as evidence for the poor 
efficiency of the energy cascade in quasiclassical QT due 
to the ‘bottleneck’ between the classical and quantum 
lengthscales m- However, recent experiments in a ro¬ 
tating container revealed vanishing traction by the con¬ 
tainer walls on turbulent superfluid '^He at low temper¬ 
atures when a < 10“^, resulting in a long-lived rotating 
state m- This cast doubt on the interpretation of ^He 
spin-down turbulence as being HIQT |26j and pushed for 
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new experiments with truly HIQT. 

Thus, the goal of this work was to measure and com¬ 
pare the decay rates of different types of turbulent flow, 
including those generated by a towed grid and impul¬ 
sive spin-down, in a broad range of temperatures. To 
determine the value of Cj one has to know both C and 
f in Q. With our technique of free decay, the injected 
energy flux, —E, is controlled by the size of the largest 
energy-containing eddy and its lifetime. In fact, 0 with 
a meaningful Q can only be applied for homogeneous tur¬ 
bulence while, for bound inhomogeneous flows, only an 
integral rate of energy dissipation can be measured to¬ 
gether with some averaged value of vortex line density. 
We will hence assume that 0 relates average E and L 
through some integral C,. 

The energy per unit mass of helium in the energy- 
containing eddies with velocity amplitude u \s E = 
where ^ ^ 1/2. Their length scale A is limited by the 
container size d, X = /3d, where /3 ~ 1. We assume that, 
as in classical turbulence, this energy is released within 
the lifetime r of order the turn-over time ^ i. e. 

T = 0Xu~^, where d ^ 1. In the quasi-steady regime, the 


energy flux fed into the cascade is hence 

—E = Et ^ or 

— 2^uu = 9~^ d~^. 

(2) 

Its solution at late time t is 

£(f) =4C6»2/32d2t-2, 

(3) 

T(t) = t/2. 

(4) 

After plugging ([^ into Q, we arrive at 

C{t) ~ Ad{Kt) , 

(5) 


where A = ~ 1. This is the C oc t 

free decay that was observed in many experiments [a EH 
Ea EH] and numerical simulations ES] . 

Our experiments were conducted in ultra-pure [301131 
^He at pressure 0.1 bar filling the volume shown in Fig.jl 
a 90 ° section of an earthed annular channel with an inner 
wall radius of curvature equal to 2.75 cm and of rectan¬ 
gular cross-section with sides dh = 1.8 cm (horizontal) 
and dv = 1.7cm(vertical). A brass grid (1.5cmx 1.5cm) 
could be electromagnetically driven at a constant velocity 
from one end of the channel to the other. The operating 
principle of the device is described elsewhere |32| while 
the technique of measuring the density and polarization 
of vortex tangles using negative ions is detailed in section 
1 of Supplemental Material (SM) [53] . 

We investigated the decay of turbulence generated by 
three different methods: towing a grid at velocity Vg ^ 
10cms“^ through the channel [38] . impulsive spin-down 
from uniform rotation at angular velocity fl ~ lrads“^ 
to rest, and injection of electric current for long periods of 
time (~ 1 nA through voltage ^ 100 V for ~ 100 s). Each 



FIG. 1. Experimental setup [32]. The front and bottom walls 
of the channel are not shown. Blue circles depict charged 
vortex rings (CVRs, not to scale), used to probe QT. CVRs 
propagate from the injector (not shown) at the front wall to 
collector inside a hole (shown by red circle) in the back wall. 
The assembly could be rotated about the vertical axis (f2). 
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FIG. 2. Decay of vortex line density jC{t) for turbu¬ 
lences, generated by different means: ‘-|-Spin-down’ - spin- 
down from Q =4-1.5 rad s“^, ‘-Spin-down’ - spin-down from 
fl=-1.5rads“^, ‘Grid 15cms“^’ - grid with rUg = 3mm and 
Vg = 15cms“^, ‘Grid 9cms“^’ - grid with rUg = 3mm and 
Vg = 9cms“^, ‘Ion jet’ - injection of negative ions at a cur¬ 
rent of 700 pA lasting for 100 s. Dashed line shows the 
dependence. T = 80 mK. 


resulted in well-developed quasi-classical turbulence in a 
wide range of length scales (the length scales and cor¬ 
responding effective Reynolds numbers are tabulated in 
section 2 of [33]). After generation, C{t) was probed 
with a pulse of ions after a delay t. Each realisation 
was probed only once to avoid distortion of the turbu¬ 
lent flow by the probing pulses. For each method, we 
forced QT sufficiently hard, that the late-time decay was 
the same, independent of the intensity of forcing (e. g. 
if Vg ^ 5cms“^). In the experiments with grid turbu¬ 
lence the values C{t) at late times did not depend on 
how many times in succession (1, 2, 3 or 10) the grid was 
towed through the channel, nor did it depend on the grid 
mesh sizes rUg used (0.75 mm and 3 mm). In the exper- 
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FIG. 3. The values of the fitting parameter A = 

vs. temperature (values of the mutual fric¬ 
tion parameter a{T) are shown at the top). Measurements 
by Stalp et al. [35] are shown by asterisks for comparison. 

iments with rotation the grid was parked at one end of 
the channel. 

For all temperatures and all methods of turbulence 
generation, after a method-specific transient process of 
duration < 10 s, the decays of vortex line density followed 
C oc as shown in Fig. We fitted them to ([^ [39] 

for time t between 30 s and 200 s, and the resulting values 
of A{T) (using d = (dh + d„)/2 = 1.75cm) are plotted 
in Fig.[^ We also compare these with the experimen¬ 
tal values of A{T) for grid turbulence (square channel, 
d= 1.27cm) |8l|35| and spin-down turbulence (cubic cell, 
d = 4.5 cm |55| and rectangular cell, d = 1.27cm p5]). 
One can see that at temperatures above 0.8 K, corre¬ 
sponding to a > 10“^, the values of A{T) for all methods 
of turbulence generation in our container agree with each 
other and also, within their scatter, with previous exper¬ 
iments. However, at T < 0.8 K A{T) approaches either 
of two zero-temperature limits: A(0) « 11 for both the 
ion-jet and grid-generated turbulence, while H(0) « 23 
for the spin-down turbulence. We would thus conclude 
that at r > 0.8 K the late-time turbulence is the same 
whatever the initial flow, i. e. approximatelly isotropic 
and homogeneous. This implies that the leftovers of the 
initial flow pattern (say, rotation following spin-down) 
disappear within less than 30 s. But at lower temper¬ 
atures, the spin-down turbulence is different from that 
for other methods at all times; this might be explained 
by our observation that the memory of initial rotation 
is retained during the late-time decay [33] - presumably 
in the form of a vortex tangle, rotating at angular veloc¬ 
ity ~ 0.1 rad s“^ near the vertical axis of the cell, that 
preserves some of the initial angular momentum. 

During the transient following the spin-down of a rect¬ 
angular cell, much of the fluid’s initial angular momen¬ 
tum is transfered to the walls through pressure fluctu¬ 
ations from large eddies, eventually creating turbulence 
with a broad distribution of length scales. At late times, 


when, as we suppose, the remains of that angular mo¬ 
mentum survive only near the axis, these pressure fluctu¬ 
ations at walls (‘pressure drag’) become inefficient, and 
only the traction at the walls (‘frictional drag’) exerts 
torque. If at T < 0.8 K this traction becomes too small 
to reduce the remaining angular momentum within the 
decay time, the effective boundary conditions (BC) be¬ 
come of the ‘slip’ type. Let us discuss two different origins 
of traction: the viscosity of the normal component and 
vortex pinning. 

For laminar flow, the relaxation time for coupling 
between the superfluid and a stationary normal com¬ 
ponent is ^ [q;(T)k£]“^. With decreasing tempera¬ 
ture, it rapidly increases and should be compared to 
the lifetime of energy-containing eddies Q: the cross¬ 
over from the limit of coupled to uncoupled components 
would thus be expected at a ~ 2[KtC{t)]~^. For typical 
C{t) ~ 10^cm“^ at t ^ 20s (as in Fig. § , this corre¬ 
sponds to a{T) ~ 10“^, i. e. to T ^ 1.1 K. However, in 
a turbulent state, the locally-enhanced density of vortex 
lines near walls might enhance the mutual friction force, 
hence allowing the cross-over to occur at smaller values of 
ct(T). Furthermore, as the mechanical forcing is expected 
to affect the large-scale superfluid and normal flow in a 
similar manner, these flows could be generated nearly 
fully-coupled from the outset; this may further ease the 
condition for coupling and allow the cross-over to ‘slip’ 
BC to occur at a lower temperature. Note that rota¬ 
tion of superfluid ^He was also found to decouple from 
container walls when a 10“ 

With numerous vortex lines terminated at the con¬ 
tainer walls, a tangential flow experiences an effective 
friction due to the pinning of these lines [ 12 ] • This force 
depends on the roughness of the surface, density of vor¬ 
tex lines as well as the lines’ dynamics - such as the 
frequency of reconnections (that facilitates effective de- 
pinning of lines) and tension in the presence of developed 
Kelvin waves. We can give a conservative estimate of the 
upper limit on this force per unit area, Fp < fC, by as¬ 
suming that all lines are strongly pinned |43j and pull in 
the direction of tangential flow with force equal to their 
line energy, / = ^ In « 1.5 pN, where Uq ~ 1A is 
the radius of the vortex core and p is the density of he¬ 
lium. Such a force would remove the angular momentum 
in a cell of square cross-section with side dh = 1.8 cm, 
initially rotating at Hq = 1.5rads“^, in ~ which 

is ^ 20 s for C{t) ~ 10"^ cm“^ at t = 20 s. While this re¬ 
laxation time is indeed comparable with the decay time, 
the force in a realistic weakly-polarized tangle should be 
much weaker. Furthermore, reconnections of pinned vor¬ 
tex lines can play an important role in promoting their 
creep from one pinning site to another |44| : this effective 
reduction of the friction force is believed to be facilitated 
by the enhanced amplitude of Kelvin waves on the scale 
of wall roughness - which are expected to rapidly grow 
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in size when a < 10“^, i. e. damping due to mutual fric¬ 
tion becomes negligible [16]. Lastly, because of frequent 
reconnections, the torque cannot extend much beyond 
one mean inter-vortex distance £q. Hence, only a vanish¬ 
ingly small shear stress can be sustained by the tangle, 
and it will be impossible to exert sufficient torque on the 
rotating core far from the container walls. 

It is thus not surprising that at T < 0.8 K, the de¬ 
coupling of the superfluid component from the container 
at large length scales and time scales of order the decay 
time results in long-lived rotation far from walls. As in 
the classical case [45], this residual rotation should slow 
down the cascade of energy to smaller eddies and thus 
increase the value of 9. Hence, according to ([^, this can 
explain the fact that, at the same decay time t, the vor¬ 
tex line density C{t) is higher for spin-down turbulence 
than for grid turbulence. It is also comforting to see that 
spin-down turbulence in containers with three different 
d all returned similar zero-temperature A = 
in Fig.(different blue symbols), as predicted by (§. 

Let us now discuss the possible effect of BC on the 
decay rate of grid turbulence. Far from walls, the dy¬ 
namics of the superfluid eddies (whether coupled to the 
low-viscosity normal component at T ^ 1 K or decou¬ 
pled from the vanishing normal component at T ^ 1K) 
at classical length scales is expected to be identical [21] . 
However, this is not the case for the energy-containing 
eddies in a container, because they are affected by walls. 
No-slip BC would speed-up the breakdown of eddies 
through the diffusion of vorticity via eddy viscosity and 
thus decrease the parameter 9 relative to its bulk value 
for eddies of the same size; while slip BC might actually 
increase the value of 9. The effective size of the largest 
eddies in a container might also be greater for slip BC, 
which will be reflected in a larger value of f3. Either 
effect could thus explain an increase of the parameter 
A oc if BC becomes of slip type below 0.8 K - 

even if the effective kinematic viscosity (£(T) stayed the 
same. 

As the values of the parameters /3 and 9 for a con¬ 
tainer of particular shape and BC are unknown, it is 
impossible to determine the accurate value of C, from 
A. Stalp et al. [5] introduced an approach, in which 
they assumed that the energy spectrum in the space 
of wavenumbers k is meaningful and equal to the Kol¬ 
mogorov spectrum Ek = (with C ~ 1.5) all 

the way down to the the cut-off wavenumber ki ^ d~^. 
In section 3 of SM [33] we show that these assumptions 
are unrealistic, and one hence cannot expect accurate 
values of C, from this approach. Yet, we quote its re¬ 
sult for T = {): for slip BC (for which ki = n/d), the 
value A(0) = 11 for grid turbulence (Fig. 4) would cor¬ 
respond to C(0) « 0.08. This agrees well with values 
C(0) = 0.08-0.09 measured experimentally [2^ ]3S] and 
C(0) = 0.06-0.10 calculated numerically [511133 for Vi- 
nen QT, in which classical degrees of freedom are not 


excited. It seems the same bulk parameter C{T) charac¬ 
terizes the efficiency of quantum cascades in HIQT for 
different spectra, thus suggesting that there is no bottle¬ 
neck between the classical and quantum cascades 

To conclude, by towing a grid through superffuid he¬ 
lium in the zero-temperature limit, we have produced 
the best-yet realization of quasiclassical HIQT filling 
a container, and measured its decay rate. The low- 
temperature decay of HIQT follows the law C oc 
observed for all quasi-classical QT, but its decay is 
markedly faster than that of the turbulence generated by 
an impulsive spin-down to rest. The latter may be due 
to the change of the effective BC from no-slip to slip be¬ 
cause of the loss of traction at the container walls below 
0.8 K. As a result, the spin-down flow maintains rotation, 
which is responsible for the slowing-down of the decay of 
turbulence. 
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